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Abstract 

The theoretical work of V.B. Braginsky predicted that radiation pressure can couple the mechan- 
ical, mirror-eigenmodes of a Fabry-Perot resonator to it's optical modes, leading to a parametric 
oscillation instability. This regime is characterized by regenerative mechanical oscillation of the 
mechanical mirror eigenmodes. We have recently observed the excitation of mechanical modes in 
an ultra-high-Q optical microcavity. Here, we present a detailed experimental analysis of this effect 
and demonstrate that radiation pressure is the excitation mechanism of the observed mechanical 
oscillations. 
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The work of V.B. Braginskyjl], 0] predicted that pressure induced by circulating radi- 
ation in a Fabry-Perot resonator can couple the optical modes to the mechanical mirror 
eigenmodes. The coupling can lead to a parametric oscillation instability, characterized by 
regenerative oscillation of the mechanical mirror eigen-modes. This mechanism has been 
studied theoretically for its possible role in setting a detection sensitivity limit in the laser 
interferometer gravitational wave observatory (LIGO) 0,Q|, but has so far not been observed 
experimentally. 

Recently, we have observed a nonlinear mechanism^ in ultra-high-Q toroid 
microcavitiesjy] that is distinct from other nonlinear mechanisms already observed in these 
structures The geometry for observation of this nonlinearity is a standard one in which 
a wave (here referred to as the pump) is coupled from a waveguide to a microcavity mode. 
The nonlinearity manifests itself as oscillations in the pump power transmitted past the 
micro-cavity. These oscillations are observed to occur at a distinct threshold pump power 
level and have spectral components at characteristic frequencies. Numerical modeling and 
spectral analysis reported in Ref. [^J revealed that the observed oscillations are due to regen- 
erative oscillation of certain mechanical eigenmodes of the toroid microcavity. In this letter, 
we demonstrate that the observed mechanical oscillations are caused by radiation pressure, 
and specifically rule out another mechanism (thermal effects j9] ). As such, this work con- 
firms the first observation of radiation-pressure-induced parametric oscillation instability as 
predicted by Braginsky. 

The theoretical treatment of Braginsky [J considered mechanical oscillations of Fabry- 
Perot mirror eigenmodes which leads to Stokes and anti-Stokes sidebands (at frequencies, 
u ±uj m where uj is the optical and u m the mechanical frequency). It was shown jl] that if the 
Stokes field coincides with an adjacent optical cavity mode the phenomenon of parametric 
oscillation instability can occur. In contrast to the Braginsky theory, we observed mechanical 
oscillations of several mechanical modes (above a certain threshold) when the mechanical 
resonance frequencies (u m ) produce Stokes and anti-Stokes fields that fall within the same 
cavity resonance (i.e. u m < ^ )[5]. For Q-factors in the range of 10 6 -10 8 this corresponds to 
frequencies in the range of ca. 1-100 MHz which coincides with the range of the first three 
fundamental mechanical modes of the toroid microcavities employed in this work. Fig. la 
shows the first, three mechanical modes of a toroid microcavity and Fig. lb their frequency 
dependence on cavity length. Note that the mechanical motion causes modulation of the 
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optical pathlength of the toroid cavity modes, causing the excitation of optical sidebands. 
These fields appear in the cavity transmission spectrum, as shown in Fig.lc. 

To account for this scenario we have extended the coupled mode analysis of Braginsky 
to the present case of Stokes and anti-Stokes frequency pairs falling within the same cavity 
resonance (For simplicity only one pair is considered here). In addition optical coupling 
effects associated with the waveguide-resonator junction are, by necessity, included in the 
analysis. Using the rotating wave and the slowly varying envelope approximation for all 
field amplitudes, the mutual coupling of the pump (a p ), Stokes (as) anti-Stokes, (a^s) and 
mechanical mode (x m ), can be described by the following coupled-mode equations: 

dx m 1 iK om . . . . 

~dT = ~2^ Xm + 2^ n C(T) ^ aAS + ^ (1) 

dctp cip . %J^fYiQ / * \ 

— = -— + iAua p + — (x m a A s + x m a s ) + us 

dt 2t ^/m^f f uj m 

da% 1 „ , . , , %K n 



— a s + % (Ad - u m ) a s — ; x m a p 



daAs 1 . / » \ iKmo 

— — = -—a As + i (Au + LO m ) a A s + — x m a p 

dt 2r ^/rn^JJUm 

In these equations, the optical pump is detuned from the cavity-mode line-center by 
Auj = uj — ujq. The Stokes and anti-Stokes frequencies lie within the resonance band- 
width of the pump mode, and, correspondingly, are detuned by Au A s = (Ao; + u m ) and 
Aa;^ = (Auj — u) m ) . The first equation describes the mechanical eigenmode where \x m \ 2 is 
normalized to mechanical energy, i.e. E = J e^aidV, (where and cr, are the diagonal 
components of the strain and stress tensor) which decays with the lifetime T m (Q m = c<j m r m ). 
Correspondingly, |a p | 2 is the energy in the pump mode, |s| 2 is the launched pump power 
in the waveguide. The total lifetime of the optical modes is given by - = — H — — , where 
the external lifetime (r ex .) describes coupling of the microcavity mode to the waveguide via 
k = i\ — and K = — is the normalized coupling constant. C(Y) is a correction factor 
[1..2] due to reduction of circulating power in the presence of modal coupling[10]. K mo = ^ 
describes the mechanically-induced displacement of the optical cavity resonant frequency 
and contains, in general, a contribution from direct spatial change as well as refractive index 
changes (stress-optical effect) 11 1 3|. The effective coupling of optical mode to the me- 
chanical mode is governed by K om = R ^ in the case of radiation pressure^]. The effective 
mass m e ff appearing in Eqn. (1) is calculated numerically, by evaluating the total mechani- 
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cal energy E m in the mechanical mode and the corresponding harmonic radial displacement 
(amplitude r) of the toroid periphery wherein the optical mode circulates (compare fig. 
4b)jl3|. Solving the coupled mode equations in steady-state the pump-power threshold for 
onset of mechanical oscillations is given by: 



w m \„2 1 1 1 1 + K + 2iAcuT 
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Pthresh — — R rn e ff — ■ — (2) 



1 r 1 1 i-i 



2r 1 + At 2 Auj 2 as 1 + 4r 2 Ao;| J 
Careful inspection of the last term of the threshold equation shows that mechanical gain is 
possible (i.e., positive threshold power) for u > uj (i.e. Auj > 0). For Auj < 0, the mechanical 
mode is damped. The need to overcome mechanical loss leads to the — dependence, 
while the dependence of radiation pressure upon circulating optical power leads to the ~ 
dependence as well as the presence of a weighting factor describing the effect of waveguide 
coupling K = ^ and pump detuning Auj. The optical-Q scaling dependences fall into 
two regimes. The first occurs when uj m < - . In this regime the mechanical oscillation 
threshold exhibits an inverse cubic dependence on optical Q (P oc ^— (7^) ) • 111 contrast, 
for uj m > - (herein called the high-frequency (HF) regime), the rapid 1/Qq dependence 
is reduced because the Stokes field build-up is less-and-less effective in creating radiation 
pressure. In this regime, minimum threshold can be shown to occur over-coupled (i.e., 
K > 1), where again the condition u m < - is met (i.e. the mechanical oscillation frequency 



is again less than t 



re "loaded" cavity bandwidth), which causes the minimum threshold 



(^' e '' dKdAu — Q to approach an asymptotic value. The transition to the HF regime, 
under conditions of optimum threshold, occurs for an optical Q given by Qq F ~ ^ e -. 

To confirm these theoretical predictions the threshold dependence (as given by Eq. 2) on 
both optical and mechanical Q-factor have been measured. The data presented are taken 
using a single microtoroid device. Coupling to the resonator was performed using a fiber- 
optic taper coupler (see inset fig.2). The micro-toroid under consideration had principal, 
pillar and minor diameters of 72, 36 and 6.8/im, respectively, and possessed mechanical 
resonances frequencies at 4.4, 25.6 and 49.8 MHz for the first three mechanical modes (n = 
1, 2, 3 and m = 0). The optical pump wavelength was ~1550 nm and mechanical oscillation 
instability was observed bv detecting the characteristic oscillations in the transmitted pump 
jower (compare Fig. lc) |5j using an electrical spectrum analyzer (ESA) as described in Ref 
Optimization of coupling (K) was performed by adjustment of the gap between the fiber 
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taper and the microtoroid as described in refs. [3, Ql- To measure the dependence of the 
oscillation threshold on Q m , a silica microprobe was brought into contact with the interior 
(disk region) of the toroid structure. Variation of the probe contact force thereby modified 
mechanical Q while leaving the optical Q unaffected. The microprobe, which was made from 
an optical fiber, had a tip diameter of ~2 fim and can be seen in the inset of Fig. 2. In 
the absence of probe contact Q m was measured to be ~ 5000 for the n = 1 mode, and upon 
progressive increase in tip pressure could be continuously decreased to Q m ~ 50. Below 
threshold, the thermal displacement of the mechanical eigenmodes (the temperature being 
300K) provides sufficient modulation to be optically detectable, causing the appearance of 
Lorenzian peaks in the cavity transmission spectrum. Q m was then determined by fitting 
the transmission spectrum with a Lorenzian, as shown in the inset of Fig. 2. For each Q m , 
the minimum threshold was measured for the n = 1 flexural mode as shown in Fig. 2. The 
solid line in the main panel shows that the data exhibit the 1/Q m dependence in agreement 
with Eqn. and Ref. Q. 

We next measured the threshold dependence on the optical Q factor as shown in Fig. 3 for 
both the n = 1 (main panel) and the n = 3 (inset) mechanical modes. The optical Q factor 
was adjusted by exciting different radial and transverse optical modes. For lower optical 
Q, wherein the acoustical oscillation frequency falls within the cavity bandwidth, the rapid 
1/Q 3 dependence is observed for n = 1 as predicted. For higher optical Q, as theoretically 
predicted a transition into the HF regime occurs at Q£ F » 10 7 . This point agrees well with 
the theoretical prediction (—). It is important to note that these observations rule out 

B v U! m ' 1 I 1 

as origin of the observed oscillations 151 ] . In Fig. 3, the solid line is the 

minimum threshold i.e. ( qkoalu = ^) as given by equation (2). With the exception of the 

effective mass, m e ff, all parameters used to create this plot were experimentally measured 

parameters (i.e., C(T), R,Q m ,uj m ,Q ,uj ). The effective mass m e ff was inferred to be rn^jf 

= 3.3 x 10~ 8 kg. 

The inset of Fig. 3 shows the measured threshold versus Q for the n = 3 mode. The 
n = 3 mode threshold dependence shows that this mode is already well into the HF regime, 
exhibiting the theoretically predicted asymptotic behavior of the minimum threshold. This 
fact is consistent with the observed resonance frequency, 49 MHz, for the n = 3 mode 
which predicts that the HF regime occurs for optical Q factors in excess of 10 7 {Qq F = 
— = 3.8 x 10 6 ). Comparison with the n = 1 mode data shows that oscillation on the 
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n = 3 mode is preferred for lower optical Qs. Indeed, preference to the n = 3 mode was 
possible by loading the microcavity into the over-coupled regime, in agreement with theory. 
The solid curve in the inset gives the single-parameter fit to the n = 3 data which yields 
n%ff = 5 x 10~ u kg, which is a factor of 660 lower than the mass of the n — 1 mode. 

As a further test of the validity of the theoretical model, the experimental effective mass 
values are compared with the theoretical prediction based on finite element modeling. For the 
n = 3 mode, the predicted mass associated with the radial motion was rrifh = 5 x 10~ n kg, 
which is in very good agreement with the experimental fit. However, for the n = 1, 2 modes, 
the calculated effective mass is a strong function of the offset of the toroidal ring with respect 
to the equatorial plane of the disk [1 61] . To both validate and quantify this offset, a cross 
section of the toroid microcavity used in this study was obtained with focused ion beam 
slicing. SEM imaging (cf. Fig. 4 panel A) reveals the presence of the above-postulated 
equatorial offset which amounts to A = 1.3 fim . Incorporation of this offset into the 
numerical mass calculation yields rn3f = 2.6 x 10~ 8 /cg and ru^ff = 2 x 10~ 9 kg. This value 
agrees very well with the experimental values from above. Finally, the numerical model 
also explains why the n = 2 mode is only observed sub-threshold in the experiments. The 
low mechanical Q value ( ~ 200) in conjunction with its high effective mass and frequency, 
predicts threshold powers > 2 mW, which are higher than pump powers available in the 
experiments. 

In summary, presented is both an experimental and theoretical analysis of radiation 
pressure induced parametric oscillation instability, as predicted by Braginsky. Excellent 
agreement of the threshold functional dependence on optical Q is obtained, providing a 
confirmation that radiation pressure is the excitation mechanism of the observed oscilla- 
tions. Besides the fundamental aspects of this work, the observed coupling of mechanical 
and optical modes by radiation pressure can find applications in micro- and nanomechanical 
systems (MEMS/NEMS) 12] for ultra-high sensitivity measurements of charge jisl ]. displace- 
ment j^|, mass, force ^ | or biological entities [^J. Equally important, radiation pressure 
as observed here can be used to achieve cooling of a mechanical resonator mode. 
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FIG. 1: Panel A: Finite element modeling of the micro-mechanical resonances of a silica toroid 
microcavity (diameter of the toroid periphery is 4 /im). The radial mode and azimuthal mode 
order are denoted with n and m, (where m = corresponds to rotationally symmetric modes). 
Shown are the first three (rotationally symmetric radial modes (n = 1, 2, 3 m = 0) in cross section 
with amplitude of motion greatly exaggerated for clarity. In addition the stress field is indicated 
using color. Note that the mechanical motion modulates the cavity pathlength due to a change in 
the cavity radius, which causes the excitation of optical sidebands. Panel B: Mechanical oscillation 
frequency versus the cavity length L for the first three fundamental mechanical modes (Dots are 
experimentally measured frequencies from Ref. I]]). Panel C: Typical frequency spectrum of the 
optical cavity transmission, revealing the presence of several sidebands which correspond to the 
first three mechanical eigenmodes (uj m ), as well as harmonics of oo m . Here, only the n = 2 mode is 
above threshold, whereas the remaining mechanical modes are subthreshold (and observable due 
to their room temperature thermal displacement noise). 
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FIG. 2: Main panel: The oscillation threshold (fiWatts) versus the mechanical quality factor 
of the n = 1 mode. The solid line is the theoretical prediction based on an inverse Q relation 
(P oc l/Qm)- Right inset: Side-view, optical micrograph of the experimental setup, consisting of a 
silica microprobe in contact with a taper-fiber-coupled micro-toroid of 72-/im-principal diameter. 
Left inset: Spectrum of the optical cavity transmission exhibiting the thermal displacement noise 
of the n = 1 mechanical mode. Solid line: The Lorenzian fit to infer the value of Q m . 
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FIG. 3: Main panel: The measured mechanical oscillation threshold (fiWatts) plotted versus the 
optical Q factor for the fundamental flexural mode (n = l,u m = 4.4 MHz,Qm = 3500). The 
solid line is a one-parameter theoretical fit obtained from the minimum threshold equation by first 
performing a minimization with respect to coupling and pump wavelength detuning, and then 
fitting by adjustment of the effective mass ( w^// = 6.6 x 10 -8 kg). Inset: The measured threshold 
for the 3 rd order mode (n = 3, v m = 49 MHz, Q = 2500) plotted versus optical Q. The solid line 
gives again the theoretical prediction with Tu^h = 1.1 x 10~ 10 kg. The n = 1 mode data from the 
main panel is superimposed for comparison. 
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FIG. 4: Panels A: Scanning-electron micrograph of the microcavity cross section achieved using 
focused-ion beam (FIB) preparation. The image reveals the presence of an offset in the toroid 
with respect to the 2-/j,m thick silica support disk (offset of ca 1.3 /im). Panel B: Finite element 
modeling of the fundamental optical mode. The presence of an offset provides a moment-arm that 
effectively enhances coupling of radiation pressure to the n = 1 transverse motion. 
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